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I Abstract 

We give necessary and sufficient conditions on a pair of positive radial functions V and W on a ball 
, B of radius R in R" , n > 1, so that the following inequalities hold for all u £ Co°(B): 

1 / s U(a;)|Vu| 2 ds > J B W(x)u 2 dx, 

-d ■ 

^2 , and 

/ B V»|Au| 2 d:c > J B W(x)\Vu\ 2 dx + (n- 1) / B (^ - ^)|Vu| 2 dz. 

This characterization makes a very useful connection between Hardy-type inequalities and the oscillatory 
behaviour of certain ordinary differential equations, and helps in the identification of a large number of 
, such couples (V, W) - that we call Bessel pairs - as well as the best constants in the corresponding 

^ ■ inequalities. This allows us to improve, extend, and unify many results -old and new- about Hardy and 

Hardy-Rellich type inequalities, such as those obtained by Caffarelli-Kohn-Nirenberg }9_, Brezis- Vazquez 
■ [S], Wang-Willem [27], Adimurthi-Chaudhuri-Ramaswamy [T], Filippas-Tertikas |13j . Adimurthi-Grossi 

ON ' -Santra [2], Tertikas-Zographopoulos [24] . and Blanchet-Bonforte-Dolbeault-Grillo-Vasquez [?]. 

7— I ' 

On ■ 

O ' 1 Introduction 

: 

Ever since Brezis- Vazquez [8] showed that Hardy's inequality can be improved once restricted to a smooth 
^ , bounded domain f2 in R™, there was a flurry of activity about possible improvements of the following type: 

Ifn>3 then f Q | Vu| 2 dx - (^) 2 / n j^dx > / n V{x)\u\ 2 dx for all u £ H^(Q), (1) 

as well as its fourth order counterpart 

If n > 5 then / n \Au\ 2 dx - "'t^ 4 ) 2 J Q ^dx > J Q W(x)u 2 dx for u £ H 2 (fl) n ff^(fi), (2) 

where F, are certain explicit radially symmetric potentials of order lower than \ (for V) and (for W) . 
In this paper, we provide an approach that completes, simplifies and improves most related results to-date 
regarding the Laplacian on Euclidean space as well as its powers. We also establish new inequalities some 
of which cover critical dimensions such as n = 2 for inequality and n = 4 for ([2]). 

We start - in section 2 - by giving necessary and sufficient conditions on positive radial functions V and 
W on a ball B in R n , so that the following inequality holds for some c > 0: 

J B V{x)\Vu\ 2 dx >cf B W{x)u 2 dx for all u £ Cfi°(B). (3) 
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Assuming that the ball B has radius R and that J Q r „-l v ^ dr — +00, the condition is simply that the 
ordinary differential equation 

(Bv.cw) y"(r) + (^ + (r) + ^»(r) = 

has a positive solution on the interval (0, R). We shall call such a couple (V, W) a Bessel pair on (0, R). The 
weight of such a pair is then defined as 

fl(V, W; R) = sup {c; {B v ^ cW ) has a positive solution on (0, i?)}. (4) 

This characterization makes an important connection between Hardy-type inequalities and the oscillatory 
behaviour of the above equations. For example, by using recent results on ordinary differential equations, 
we can then infer that an integral condition on V, W of the form 

lunmpr 2 ^V(r)W(r)( j* *L— f < 1 (5) 



r— >0 



V(t) 



is sufficient (and "almost necessary") for (V, W) to be a Bessel pair on a ball of sufficiently small radius p. 

Applied in particular, to a pair (V, p-V) where the function is assumed to decrease to —A on (0, R), 

we obtain the following extension of Hardy's inequality: If A < n — 2, then 

J B V(x)\Vu\ 2 dx > i 11 ^ 2 -) 2 J B V(x)^dx for all ue C* °°(B) (6) 

and ( "~2~ 2 ) 2 is the best constant. The case where V(x) = 1 is obviously the classical Hardy inequality and 
when V(x) = |x|~ 2a for — 00 < a < £i 2 l2 -, this is a particular case of the Caffarelli-Kohn-Nirenberg inequality. 
One can however apply the above criterium to obtain new inequalities such as the following: For a, b > 

• If a(3 > and m < 2=2 then for all u e C$°(R n ) 

(a + feM a )V ,2 , . ( n-2m-2 f (a + b\x\^ 

lx{2m \Vu\ dx > ( ) jf^ |a;|2m+2 u dx, (7) 

and ( "~ ^ ~ ) 2 i s the best constant in the inequality. 

• If a/3 < and 2m - a/3 < n - 2, then for all u € C£°(R") 
{a + b\x\ a f 2 n- 2m + a/3-2 2 /" (a + b\x\ a f 2 



and ( "~ 2m +"^~ 2 )2 j g Des j. cons |; a nt in the inequality. 
We can also extend some of the recent results of Blanchet-Bonforte-Dolbeault-Grillo-Vasquez [4]. 
• If a/3 < and -a/3 < n - 2, then for all u e C$°(M. n ) 

{a + b\x\ a f\X7u\ 2 dx>bi{ U ~ a f~ 2 ) 2 [ (a + b\x\ a f-^u 2 dx, (9) 



and b& ( " a ^ 2 ) 2 is the best constant in the inequality. 
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• If a/3 > 0, and n > 2, then there exists a constant C > such that for all u € C^(R™) 

(a + blxn^Vwpcfe > C [ (a + b\x\ a f~iu 2 dx. (10) 

n-2\2 ^ n S U- fn+af3-2\2 



Moreover, {^) 2 < C < 6» ( , 
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On the other hand, by considering the pair 

V(x) = \x\~ 2a and W a , c (x) = { n^l f\ x \-^-2 + c\x\- 2a W(x) 
we get the following improvement of the Caffarelli-Kohn-Nirenberg inequalities: 

f \x\~ 2a \\7u\ 2 dx- ( n ~ 2 °~~ 2 ) 2 / \x\- 2a ~ 2 u 2 dx>c [ \x\- 2a W{x)u 2 dx for all we C °°(B) (11) 

if and only if the following ODE 

(BcwO v" + \v' + cW ( r )y = 

has a positive solution on (0, R). Such a function W will be called a Bessel potential on (0, R). This type of 
characterization was established recently by the authors |15j in the case where a = 0, yielding in particular 
the recent improvements of Hardy's inequalities (on bounded domains) established by Brezis- Vazquez [8], 
Adimurthi et al. pQ, and Filippas-Tertikas [13]. Our results here include in addition those proved by Wang- 
Willem [57] in the case where a < and W[r) — r2 ^Ry i but also cover the previously unknown limiting 

case corresponding to a = as well as the critical dimension n = 2. 

More importantly, we establish here that Bessel pairs lead to a myriad of optimal Hardy- Rellich inequalities 
of arbitrary high order, therefore extending and completing a series of new results by Adimurthi et al. [2J, 
Tcrtikas-Zographopoulos [24j and others. They are mostly based on the following theorem which summarizes 
the main thrust of this paper. 

Theorem 1.1 Let V and W be positive radial C 1 -functions on B\{0}, where B is a ball centered at zero 
with radius R in R" (n > 1) such that rn -i v ^ dr = +oo and j^' r n ~~ 1 V(r)dr < +oo. The following 
statements are then equivalent: 

1. (V, W) is a Bessel pair on (0, R) and (3(V, W; R) > 1. 

2. J B V{x)\Vu\ 2 dx > f B W{x)u 2 dx for all u G Cg°(B). 

3. If hnv^o r a V(r) = for some a < n — 2, then the above are equivalent to 

J B V{x)\Au\ 2 dx > J B W{x)\Vu\ 2 dx + (n - 1) ^(^p 1 - ^^-)\Vu\ 2 dx for all radial u € C^ r (B). 

4- If in addition, W(r) — 2V r ^ + 2Vr J r ^ — V rr (r) > on (0,R), then the above are equivalent to 

J B V{x)\Au\ 2 dx > J B W{x)\Vu\ 2 dx + (n - 1) f B (^ - ^^)\Vu\ 2 dx for all u € C^{B). 

In other words, one can obtain as many Hardy and Hardy-Rellich type inequalities as one can construct 
Bessel pairs on (0, R). The relevance of the above result stems from the fact that there are plenty of such 
pairs that are easily identifiable. Indeed, even the class of Bessel potentials -equivalently those W such that 
(l, ( 23 ^ L =) 2 2 + cW(x)j is a Bessel pair- is quite rich and contains several important potentials. Here are 
some of the most relevant properties -to be established in an appendix- of the class of C 1 Bessel potentials 
W on (0, R), that we shall denote by B(0, R). 

First, the class is a closed convex solid subset of C 1 (0, R), that is if W G B(0, R) and < V < W, then 
V e B{0, R). The "weight" of each W G B(R), that is 

[3(W; R) — sup {c > 0; (B c w) has a positive solution on (0, i?)}, (12) 

will be an important ingredient for computing the best constants in corresponding functional inequalities. 
Here are some basic examples of Bessel potentials and their corresponding weights. 

• W = is a Bessel potential on (0, R) for any R > 0. 
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• W = 1 is a Bessel potential on (0, R) for any R > 0, and R) = ^ where zq — 2.4048... is the first 
zero of the Bessel function Jo • 

• If a < 2, then there exists R a > such that W(r) = r~ a is a Bessel potential on (0,i? a ). 

,e((t-l)-Mme.) 

• For k > 1, R > and p = R(e e " ), let W k , p (r) = S^ =1 ^(n*=i W ?)" 2 where the func- 
tions Zo<?W are defined iteratively as follows: log^ x \.) = log(.) and for k > 2, log^(.) = log{log^ k ^ x \.)) . 
Wk^ is then a Bessel potential on (0, R) with (3(Wk. p ; R) = 

• For jfe > 1, i? > and p > R, define W4 ;p (r) = S* = i^i(p^I(p • • • X j-d^) X j (£) where the Ac- 
tions Xj are defined iteratively as follows: Xi(t) = (1 — log(t)) -1 and for fc > 2, X k (t) = X 1 (Xfc_ 1 (t)). 
Then again Wfc iP is a Bessel potential on (0,R) with f3{Wf., p ',R) = t- 

• More generally, if W is any positive function on M such that liminf ln(r) JT sW(s)ds > — oo, then for 
every R > 0, there exists a := a(-R) > such that W a {x) := a 2 W(ax) is a Bessel potential on (0,R). 

What is remarkable is that the class of Bessel potentials W is also the one that leads to optimal improvements 
for fourth order inequalities (in dimension n > 3) of the following type: 

J B \Au\ 2 dx-C(n)J B l -^fdx>c(W,R)J B W(x)\Wu\ 2 dx for all u E H 2 (B), (13) 

where C(3) = ||, C(4) = 3 and C(n) = ^- for n > 5. The case when W = Wk. P and n > 5 was recently 
established by Tertikas-Zographopoulos [24] . Note that W can be chosen to be any one of the examples of 
Bessel potentials listed above. Moreover, both C(n) and the weight f3(W; R) are the best constants in the 
above inequality. 

Appropriate combinations of ^ and p3[) then lead to a myriad of Hardy-Rellich inequalities in dimension 
n > 4. For example, if W is a Bessel potential on (0, R) such that the function r decreases to —A, and 
if A < n — 2, then we have for all u E C^(B^) 

IA „, 7i 2 (n-4) 2 f u 2 , ,n 2 (n-A-2) 2 N . f W(x) a , 

By applying (|14[) to the various examples of Bessel functions listed above, one improves in many ways the 
recent results of Adimurthi et al. [5] and those by Tertikas-Zographopoulos [M]- Moreover, besides covering 
the critical dimension n = 4, we also establish that the best constant is (1 + "(" g ~ 4 ) ) for all the potentials 
Wk and Wk defined above. For example we have for n > 4, 

L >- nJ ^L ^'gX wW"^ <»> 

More generally, we show that for any to < ^^j^, and any W Bessel potential on a ball Br C i? n of radius R, 
the following inequality holds for all u E C^(B^) 

IAd 2 f IVmI 2 f iVul 2 

/;B |-p>«n, m y^^p^^+^i?) y^^)^^, (i6) 

where a m . n and R) are best constants that we compute in the appendices for all to and n and for many 
Bessel potentials W. Worth noting is Corollary 13.21 where we show that inequality (fTB]) restricted to radial 
functions in C§°(Br) holds with a best constant equal to ( n+2m ) 2 , but that a n<m can however be strictly 
smaller than ( " + 2 2m ) 2 in the non-radial case. These results improve considerably Theorem 1.7, Theorem 1.8, 
and Theorem 6.4 in |2~4"1. 
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We also establish a more general version of equation (fl4|) . Assuming again that r ^ff decreases to —A on 
(0, R), and provided m < and A < n — 2m — 2, we then have for all u € Cg°(Bn), 

f ^dx > 3 f -^dx 

w (n + 2rn) 2 (n-2m-A-2) 2 , f W(x) 
+P(W;R)( [ 4 > + [ - >-) j^u 2 dx, (17) 

where again the best constants B n ^ m are computed in section 3. This completes the results in Theorem 1.6 
of [24], where the inequality is established for n > 5, < m < and the particular potential Wk, p - 
Another inequality that relates the Hessian integral to the Dirichlet energy is the following: Assuming 
-1 < m < 2=± and W is a Bessel potential on a ball B of radius R in R n , then for all u £ Cq°(B), 

[ ^- (n + 2m)2( V 2m ~ 4)2 / r^dx > 3(W;R)^ + y 2 [ ^u 2 dx 
J B \x\ 2m 16 Jb \x\ 2m + 4 - ^ V ' 4 J B |z| 2 ™+ 2 

+B(\x\ 2m ] R)\\u\\ Hh . (18) 

This improves considerably Theorem A. 2. in [2] where it is established - for m — and without best 
constants - with the potential Wi tP in dimension n > 5, and the potential W2. P when n = 4. 

Finally, we establish several higher order Rellich inequalities for integrals of the form J Br dx, improving 

in many ways several recent results in [24] . 

The approach can also be used to improve the recent results of Liskevich-Lyachova-Moroz [TH] on exterior 
domains and will be developed in a forthcoming paper. 

2 General Hardy Inequalities 

Here is the main result of this section. 

Theorem 2.1 Let V and W be positive radial C 1 -functions on Bf>\{0\, where Br is a ball centered at zero 
with radius R (0 < T < +oo) in K™ (n > 1). Assume that f£ rn -i v , r \ dr = +oo and r n ~ 1 V{r)dr < oo 
for some < a < R. Then the following two statements are equivalent: 

1. The ordinary differential equation 

(Bv.w) y"{r) + (2=1 + ^g)y'(r) + %$y(r) = 

has a positive solution on the interval (0, R) (possibly with ip(R) = 0). 

2. For allue C$°(B R ) 

(R VtW ) J Br V(x)\Vu(x)\ 2 dx > J Br W(x)u 2 dx. 

Before proceeding with the proofs, we note the following immediate but useful corollary. 

Corollary 2.2 Let V and W be positive radial ^-functions on B\{0}, where B is a ball with radius R in 
R" (n > 1) and centered at zero, such that f Q r „-l v , r \ dr = +oo and r n ~ 1 V{r)dr < oo. Then (V, W) is 
a Bessel pair on (0, R) if and only if for all u £ C^'(Br), we have 

I V{x)\X7u\ 2 dx> 0(V 7 W;R) f W(x)u 2 dx, 
Jb r Jb r 

with 3(V, W] R) being the best constant. 

For the proof of Theorem 12. 1( we shall need the following lemmas. 
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Lemma 2.3 Let Q, be a smooth bounded domain in R n with n > 1 and let ip G C ll (0,i? := sup^ggo, \ x \) ^ e a 
positive solution of the ordinary differential equation 



„ , ,n-l , F r (r). , , W(r) 



y +( + T^i )v + ttHv = o. ( 19 ) 

on (0,i?) /or some V(r),W(r) > where ^-ly^ dr = +oo and J Q r n ~ 1 V(r)dr < oo. Setting i]{x) — 
/or any u € Co°(fi) ; we £/ien /iaue £/ie following properties: 

1. ^r n ^V{r){^fdr<oo and lim r ^ r^MO^ = 0. 
«. / F(|x|)(^'(|x|))V 2 (a;)^ <oo. 
5. / n F(|x|)<p 2 (|a;|)|VVf (a;)cfe < oo. 

4. |/ n nklV(klMI^W^-VV(^|<oo. 

5. lim,.^o I ^(M)v?'(l a; l)v?(l a; l)' ( /' 2 ( a: ) a ' s l = 0, where B r <Z Q, is a ball of radius r centered at 0. 
Proof: 1) Setting x(r) = r n_1 V(r)^pj- , we have 

r-V>K(r) + * 2 (r) = I^^(^ (r)+( !Lzl + ^y (r)) = . ^(O^) < 0; o < r < -R. 

r- v(r) 

Dividing by r"~ 1 V A (r) and integrating once, we obtain 

X{r) ^[^%) dS + X{R) - (20) 

To prove that lim r ^o G(r) < oo, where G(r) :— rfs, we assume the contrary and use ([201) to write 

that 

(_ r »-iy( r ))G'(r))5 > G(r) + x(R). 

Thus, for r sufficiently small we have — r" -1 ^?-)^'^) > ^G 2 (r) and hence, (-Qj-pj)' > 2r"-*y(r) ' wn i cn 
contradicts the fact that G(r) goes to infinity as r tends to zero. 

Also in view of (|20p . we have that Xq := lim r ^o x ( r ) exists, and since lim r _>o G{r) < oo, we necessarily have 

xq = and 1) is proved. 

For assertion 2), we use 1) to see that 

V(\x\W(\x\)) 2 ^(x)dx<\\u\t [ V(\x\) { -^M^dx<oo. 

Jn 9 (fI) 

3) Note that 

< + K*)I^S < + C 2 l$H§, for all x e CI, 

where C\ = max ie n |Vm| and C 2 = max xS n \u\. Hence we have 

(c^(M) + cv(M)) 2 



f V(\x\)^(\x\)\V^\ 2 (x)dx < f V(\x\)- 



<p*(\x\) 



dx 



C 2 V(\x\)dx + f 2C 1 C^^V{\x\)dx+ [ Cl(^-fV{\x\)da 
Jn <P{\ X \) Jn ¥>(M) 

< L 1 + 2C 1 C 2 ( [ V{\x\)C^^-fdxY{ [ V{\x\)dxY+L 2 

Jn f(\ x \) Jn 

< 00, 



G 



which proves 3). 

4) now follows from 2) and 3) since 

V(\x\)\Vu\ 2 = V{\xW{\x\)ftf{x) + 2V(\x\)^'(\x\M\x\)i,(x)-^.V^(x) + V (\x\)<p 2 (\x\)\V^ 

Finally, 5) follows from 1) since 

V(\x\)<p>(\x\M\x\)^(x)ds\ < \\u\a [ V(\x\)^^-ds 

,</;, JdB r vim) 

2 >-,J<P'(r)\ 



,2 „n-lr,^k'WI 



</?(r) 



Lemma 2.4 Let V and W be positive radial C 1 -functions on a ball B\{0}, where B is a ball with radius R 
in R" (n > I) and centered at zero. Assuming 

J B (V{x)\Vu\ 2 - W{x)\u\ 2 ) dx>0 for all u € Cg°(B), 

£/ien i/iere exists a C 2 -super solution to the following linear elliptic equation 

- div(V(x)Vu) - W{x)u = 0, in B, (21) 
u > in 5\{0}, (22) 
u = in dB. (23) 

Proof: Define 

^ M ■ ; J B nx)\^\ 2 -W(xM 2 

Ai(V):=mf{ j |^| 2 ; ^eC (5\{0})}. 

By our assumption Ai(V) > 0. Let (<p n , A") be the first eigenpair for the problem 

(L — Ai(y) — Xi)ip n = on B\Br 

ip n = on 8{B\Br), 

where Lu = — div(V(x) Vu) — W(x)u, and B« is a ball of radius ^, n > 2 . The eigenfunctions can be 

chosen in such a way that ip n > on B \ Br and ip n (b) — 1, for some b £ B with < \b\ < R. 

Note that A" J, as n — > oo. Harnak's inequality yields that for any compact subset K, < C(K) 

with the later constant being independant of tp n . Also standard elliptic estimates also yields that the family 
(ip n ) have also uniformly bounded derivatives on the compact sets B — Br. 

Therefore, there exists a subsequence (<Pm 2 )i 2 °f {fn)n such that {<Pm 2 )i 2 converges to some cf2 € C 2 (B\ 
B(~)). Now consider (<Pn l2 )i 2 on B\ ^("f)- Again there exists a subsequence (<Pm 3 )l 3 01 (fni 2 )i 2 which 
converges to cp 3 E C 2 (B \ B(^)), and ips(x) = f2{x) for all x G -B \ -B(S). By repeating this argument we 
get a supersolution tp E C 2 (B\ {0}) i.e. L</? > 0, such that y> > on B \ {0}. □ 



Proof of Theorem I2.lt First we prove that 1) implies 2). Let cp € C 1 (0,i?] be a solution of such 

u(x) 

Vu\ 2 = fo/flsDjVfc) + 2^(|:r|M|z|)V>(x)^.W> + ^(|a:|)|VV| 5 



that <p(a;) > for all x € (0, i?). Define = ^(x). Then 



Hence. 



V(\x\)\Vu\ 2 > V{\xW{\x\)f^{x) + 2V(\x\)<p'(\x\M\x\)i,(x)^W(x). 
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Thus, we have 

/ V{\x\){ V '{\x\)) 2 ^ 2 {x)dx + / 2V{\x\)ip'{\x\)<p{\x\)TP{x) 1 ^ 1 .Vipdx. 

Jb Jb m 

Let B e be a ball of radius e centered at the origin. Integrate by parts to get 

^(|izi|)|Vu| 2 (fa > / V(\x\)(ip' {\x\)) 2 ip 2 (x)dx + / 2V(\x\)ip\\x\)ip(\x\U(x)-^ .Vipdx 
Jb Jb. fI 

2V(\x\W(\x\)<p{\x\)il>(x)-Z T .Vil>dx 

b\b c m 

V(\x\){<p'{\x\)) 2 ip 2 {x)dx+ [ 2V(\x\)<p'(\x\M\x\)1>(x)-£r.Vil>dx 



B\B C 



{V(\x\)<p"(\x\M\x\) + ( {n + V P (|x|)y (HMIa;!))^ 2 ^) |> d.r 

V(\x\)v'(\x\M\x\)il> 2 (x)ds 

(S\Be) 

Let e -> and use Lemma |2~51 and the fact that ip is a solution of (D VtW ) to get 

ni^l)|Vu| 2 dx > - / [nHV"(H) + ( ( "~ lMN) +K.(N)V(N)]^rfx 

r tp[\x\) 

2/ 



W(|x|)w 2 (a;)da;. 

To show that 2) implies 1), we assume that inequality (Hy,w) holds on a ball B of radius R, and then apply 
Lemma to obtain a C 2 -supersolution for the equation (|2~lj) . Now take the surface average of u, that is 



nuj. 



— - { u(x)dS = / u(ruj)d 

>r n JdB r nuj n J H=1 



V( r ) = / u ( x ) dS = / u(rcu)dLO > 0, (24) 



where oj n denotes the volume of the unit ball in R n . We may assume that the unit ball is contained in B 
(otherwise we just use a smaller ball). We clearly have 

y » {r) + ^-ly'(r) = — * f Au(x)dS. (25) 
r nLo n r JdB r 

Since u(x) is a supersolution of (|2Tj) . we have 



and therefore, 



It follows that 



div(V(\x\)Vu)ds - / VK(|a;|)wdx > 0, 



F(r) / Awc^S - K(r) / Vu.xds - W(r) / u(x)ds > 0. 

JSB r JSB r JdB r 



V(r) / AudS - V r (r)y'(r) - W{r)y{r) > 0, (26) 

JdB r 

and in view of (flM)) , we see that y satisfies the inequality 

V(r)y"(r) + (^—^M + 7 r (r))y'(r) < -W(rWr), /or < r < R, (27) 
r 

that is it is a positive supersolution for (By,w)- 

Standard results in ODE now allow us to conclude that (Bv } w) has actually a positive solution on (0,R), 
and the proof of theorem 12. II is now complete. 



2.1 Integral criteria for Bessel pairs 

In order to obtain criteria on V and W so that inequality (H^wO holds, we clearly need to investigate 
whether the ordinary differential equation (By,w) has positive solutions. For that, we rewrite (By,w) a s 

(r^V^y'y + r^W^y = 0, 

and then by setting s = - and x(s) = y(r), we see that y is a solution of (By.w) on an interval (0, 5) if and 
only if x is a positive solution for the equation 

( s -(«-3)y(iy( s ))/ + s -(n+i) W (I) x ( s ) =0 on (i )0O ). (28) 

Now recall that a solution x(s) of the equation ([28]) is said to be oscillatory if there exists a sequence {a n }^ = i 
such that a n — > +(X) and x(a n ) = 0. Otherwise we call the solution non-oscillatory. It follows from Sturm 
comparison theorem that all solutions of (j2"5|) are either all oscillatory or all non-oscillatory. Hence, the fact 
that (V, W) is a Bessel pair or not is closely related to the oscillatory behavior of the equation ([28]) . The 
following theorem is therefore a consequence of Theorem 12.11 combined with a relatively recent result of 
Sugie et al. in [22] about the oscillatory behavior of the equation (|28|) . 



Theorem 2.5 Let V and W be positive radial C 1 -functions on Br\{0}, where Br is a ball centered at 
Assume 



with radius R in W 1 (n > X). Assume f Tn ^l v , T \ dT — +oo and r™ 1 v{r)dr < oo. 



umsup r 2 ^V(r)W(r) ( J Tn J v{r) 2 < \ (29) 

f/ien (V, W) is a Bessel pair on (0,p) for some p > and consequently, inequality (Rv,w) holds for all 
u 6 C^°(B p ), where B p is a ball of radius p. 

• On i/ie o£/ier hand, if 

limMr^V(r)W(r)(J* ' —L_dr) 2 > \ (30) 

then there is no interval (0, p) on which (V, is a Bessel pair and consequently, there is no smooth 
domain f2 on which inequality (Jiy,w) holds. 

A typical Bessel pair is (|a;|~ A , |x|~ A ~ 2 ) for A < n — 2. It is also easy to see by a simple change of variables 
in the corresponding ODEs that 

W is a Bessel potential if and only if (|x| _A , |x| _A (|x| -2 + W^dxD) is a Bessel pair. (31) 

More generally, the above integral criterium allows to show the following. 

Theorem 2.6 Let V be an strictly positive C 1 -function on (0, R) such that for some A e M 

+ A > on (0, R) and lim r -^$- + A = 0. (32) 

If X < n — 2, then for any Bessel potential W on (0,R), and any c < f3(W;R), the couple (V,W\ tC ) is a 
Bessel pair, where 

W x ,c(r) = V(r)( C'^ 2 ) 2 r- 2 + cW(r)). (33) 
Moreover, 0(V, W\, c ; R) = 1 for all c < (3(W; R). 
We need the following easy lemma. 
Lemma 2.7 Assume the equation 

y" + V + V{r)y = 0, 
r 

has a positive solution on (0,R), where a > 1 and V(r) > 0. Then y is strictly decreasing on (0, R). 
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Proof: First observe that y can not have a local minimum, hence it is either increasing or decreasing on 
(0, 6), for 5 sufficiently small. Assume y is increasing. Under this assumption if y'{a) — for some a > 0, then 
y"(a) = which contradicts the fact that y is a positive solution of the above ODE. So we have ^jr < —f, 
thus, 

Therefore, x{r) — > — oo as r — > which is a contradiction. Since, y can not have a local minimum it should 
be strictly decreasing on (0, R). □ 
Proof of Theorem I2.6t Write -wrr = — - + f(r) where f(r) > on (0, R) and lim rf(r) = 0. In order to 
prove that (V(r), V(r)(( n ^^ 2 ) 2 r~ 2 + cW(r))j is a Bessel pair, we need to show that the equation 

y" + ( n ~ A ~ 1 + f(r))y' + (( n ~^~ 2 ) 2 r- 2 + cW(r))y{r) = 0, (34) 
r 2 

has a positive solution on (0,i?). But first we note that the equation 

x " + ( n ~ A ~ 1 y + (( n ~ A ~ 2 ) 2 r~ 2 + cW(r))x(r) - 0, 
r 2 

has a positive solution on (0, R), whenever c < j3(W;R). Since now /(r) > and since, by the proceeding 
lemma, x'(r) < 0, we get that x is a positive subsolution for the equation (|34|) on (0, R), and thus it has a 
positive solution of (0, R). Note that this means that f3(V, W\ jC ; R) > 1. 

For the reverse inequality, we shall use the criterium in Theorem l2.5l Indeed apply criteria (|29| to V(r) and 
W 1 (r) = C^ to get 

R i /--R i 



li mr 2(«-iV( r )VF !(r)( / TTTTT^Y = Climr 2 ("- 2 V 2 (r)( / TTTTT^) 

r^O n Jr T n ~ l V(T) ' r^O K Jr T^V (j) ' 



C( lim r("- 2 V(r) / TTTTT^) 



R , 

1 , \2 



C( lim 



r ^ (n-2)r"- 3 y(r)+r"- 2 y r (r) 

r 2(n-2)y2( r ) 

C( lim - — tttt) 2 

V -° (n _ 2)+r ^ 



(n — A - 2) 2 ' 

For (V,CV(r~ 2 + cW)) to be a Bessel pair, it is necessary that ^ n _^_ 2 yi < j, an d the proof for the best 
constant is complete. □ 
With a similar argument one can also prove the following. 

Corollary 2.8 Let V and W be positive radial ^-functions on Br\{0}, where Br is a ball centered at zero 
with radius R in W l (n > I). Assume that 

lim r X r S r ) = —A and A < n — 2. (35) 

• // lim sup r 2 '^Vf\ ^ ( "~o~ 2 ) 2 > then (V, W) is a Bessel pair on some interval (0,p), and consequently 
there exists a ball B p C R n such that inequality (Hv,w) holds for all u £ C^°(B P ). 

• On the other hand, if lim inf r 2 ^S r ) > ( "~^~ 2 ) 2 . then there is no smooth domain £1 C R n such that 

' J r ^Q v ( r > 2 

inequality (Hv.w) holds on CI. 
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2.2 New weighted Hardy inequalities 

An immediate application of Theorem 12.61 and Theorem 12.11 is the following very general Hardy inequality. 

Theorem 2.9 Let V(x) = V(|x|) be a strictly positive radial function on a smooth domain O containing 
such that R = sup x£ Q \x\. Assume that for some A £ R 

r -^- + A > on (0, R) and lim r -^- + A = 0. (36) 

1- If A < n — 2, then the following inequality holds for any Bessel potential W on (0, R): 

J n V(x)\Wu\ 2 dx > (2^2) 2 f n ¥^- u 2 dx + 0(W;R) J n V{x)W(x)u 2 dx for allueCg°(n), (37) 

and both ( "~^~ 2 ) 2 and fl(W;R) are the best constants. 
2. In particular, (i{V,r~ 2 V; R) = ( "~^~ 2 ) 2 is the best constant in the following inequality 

j n V(x)\Vu\ 2 dx> (i^ 2 ) 2 J n ^u 2 dx for allu€CZ°(tt). (38) 

Applied to Vi(r) = r- m W k Jr) and V 2 (r) = r- m W Kp {r) where W Kp {r) = E* =1 £(IIi=i log®* and 
Ty fe;p (r) = Ej =1 ^-X 2 (^)A|(^) . . . X?_ 1 (-)X?(-) are the iterated logs introduced in the introduction, and 
noting that in both cases the corresponding A is equal to 2to + 2, we get the following new Hardy inequalities. 

Corollary 2.10 Let Q be a smooth bounded domain in M. n (n > 1) and m < Then the following 
inequalities hold. 

^4^\Vu\ 2 dx > ( n ' 2 o m ' 4 ) 2 / ^Mu 2 dx (39) 

\x\ 2m 11 ~ v 2 ' J n \x\ 2m + 2 K ' 

, ^4^\Vu\ 2 dx > ( n ~ 2 o m - 4 ) 2 / P^u 2 dx. (40) 
In \x\ 2m 1 1 " V 2 ' J n \x\ 2m + 2 v ' 

Moreover, the constant ( "~ 2 "~ 4 ) 2 is £/ie &esi constant in both inequalities. 

Remark 2.11 The two following theorems deal with Hardy- type inequalities on the whole of R™. Theorem 
12.11 already yields that inequality {Hy,w) holds for all u € C^°(lR ra ) if and only if the ODE (Bv,w) has a 
positive solution on (0, oo). The latter equation is therefore non-oscillatory, which will again be a very useful 
fact for computing best constants, in view of the following criterium at infinity (Theorem 2.1 in [22j ) applied 
to the equation 

(a(r)y')' + b(r)y(r) =0, (41) 

where a(r) and b(r) are positive real valued functions. Assuming that -^fyjdr < oo for some d > 0, and 
that the following limit 

(/•oo y x 2 

/ rfr 
Jr a ( r ) 

exists. Then for the equation (|4ip equation to be non-oscillatory, it is necessary that L < i. 
Theorem 2.12 Let a,b> 0, and a,{3,m be real numbers. 



If a (3 > 0, and m < then for all u £ Cfi°(} 



<!±5g2W > { i^zi f i a+M^u'^ (4 2 ) 



' .Ln \X\ 2 " 1 + 2 

and ( "~ 2 I"~ 2 ) 2 is the best constant in the inequality 
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• If a/3 < 0, and 2m - a(3 < n - 2, then for all u € Cg°(R n ) 



and ^ "~ 2m + a / 3 ~ 2 )2 j s ^ e ^ es -j- cons i an i i n ffe inequality. 



Proof: Letting V(r) = {a+b 2 r P\ then 



V'(r) „ 6a/3r Q „ „ aa/3 

• w =-2mH -- — = -2m + a/3- 



V(r) a + br a a+br a 



Hence, in the case a, (3 > and 2m < n — 2, (|42|) follows directly from Theorem 12.91 The same holds for 
(|43|) since it also follows directly from Theorem 12.91 in the case where a < 0, (3 > and 2m — a/? < n — 2. 
For the remaining two other cases, we will use Theorem 12.11 Indeed, in this case the equation (By,w) 
becomes 

,, n— 2m— 1 6a/3r Q_1 , , c 

2/ + ( + 1 7 - V + -y = 0, (44) 

r a + br a r z 

and the best constant in inequalities (|42|) and (j43|) is the largest c such that the above equation has a positive 
solution on (0, +oo). Note that by Lemma [2~71 we have that y' < on (0, +oo). Hence, if a < and < 0, 
then the positive solution of the equation 

~ i /n— 2m — 2\2 

n — 2m — 1 



-v + - — = 

i r A 

is a positive super-solution for (j44|) and therefore the latter ODE has a positive solution on (0, +oo), from 
which we conclude that |42|) holds. To prove now that ( "~ 2 ™~ 2 ) 2 is the best constant in (|4*2"j) . we use the fact 
that if the equation (|44[) has a positive solution on (0, +oo), then the equation is necessarily non-oscillatory. 
By rewriting (|44]l as 

(r n - 2m - X (a + br a y 3 y')' + cr n - 2m - 3 (a + br a fy = 0, (45) 

and by noting that 

/>oo 

< OO, 



d r n ~ 2m - 1 (a + br a )P 



and 

lim cr 2 ("- 2m - 2 )(a + 6r Q ) 2/3 ( f°° = ^ = 

r-^oo v ; \J r r"- 2m - 1 (a + 6r«)' 3 / (n-2m-2) 2 ' 

we can use Theorem 2.1 in [32] to conclude that for equation (|4"5]l to be non-oscillatory it is necessary that 



c 1 

< 



(n- 2m- 2) 2 ~ 4 



Thus, — 2 ™ ^ is the best constant in the inequality (|4"2"| . 
A very similar argument applies in the case where a > 0, [3 < 0, and 2m < n — 2, to obtain that inequality 
flUJl holds for all u € Cg°(R") and that ( n ~ 2m + a P- 2 f i s indeed the best constant. □ 

Note that the above two inequalities can be improved on smooth bounded domains by using Theorem 12.91 
We shall now extend the recent results of Blanchet-Bonforte-Dolbeault-Grillo-Vasquez [3] and address some 
of their questions regarding best constants. 

Theorem 2.13 Let a,b > 0, and a,f3 be real numbers. 

• Ifa(3 < and -a/3 < n - 2, then for all u E C£°(R n ) 



{a + b\x\ a f\Wu\ 2 dx>b-{ £ ) 2 / (a + b\x\ a f--u 2 dx, (46) 

and j)a( "~ a ^~ 2 ) 2 ffte best constant in the inequality. 
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• If a (3 > and n > 2, then there exists a constant C > such that for all u € C^°(R n ) 

{a + b\x\ a f\Vu\ 2 dx > cj {a + b\x\ a f-^u 2 dx. (47) 

Moreover, fo§(^) 2 < C < ftl ( n + a ^- 2 f , 
Proof: Letting V(r) = (a + frr")' 3 , then we have 

V'(r) ba(3r a aa/3 
r — — — = — = a(3 — 



V(r) a + br a a + br a ' 

Inequality (|46|) and its best constant in the case when a < and /3 > 0, then follow immediately from 
Theorem 12. 91 with A = —a/3. The proof of the remaining cases will use Theorem 12. II as well as the integral 
criteria for the oscillatory behavior of solutions for ODEs of the form {By,w)- 

Assuming still that a(3 < 0, then with an argument similar to that of Theorem l2.12l above. one can show that 

-af3- 
4 r 2 



the positive solution of the equation y" + ( n+al3 1 )y' + ^ n+a f 2 ^ y = 0on (0, +oo) is a positive supersolution 



for the equation 



„ ,n-l V'(r) , b^{n + ap~2f 

y + + tttt y + — - — : — -i — y = o. 

r V(r) A(a + br a )- 



TheoremOthen yields that the inequality (J46J) holds for all u e C^°(R"). To prove now that bi ( n+( f~ 2 ) 2 
is the best constant in (|4"o]l it is enough to show that if the following equation 

(r"- x (a + br a fy')' + cr^^a + br a f~iy = (48) 

has a positive solution on (0, +oo), then c < bi ( " + ° /9 " 2 ) 2 . If now a > and < 0, then we have 

lim cr 2{n - 1] (a + br a ) 2f} -« ( [°° - dA = - 

r^oo \J r r"- 1 (a + 6r«)' 3 J hi (n + a(3 - 2) 2 ' 

Hence, by Theorem 2.1 in 22] again, the non-oscillatory aspect of the equation holds for c < b ° ("+^~ 2 ) 
which completes the proof of the first part. 

A similar argument applies in the case where a/3 > to prove that (|4"T)) holds for all u E Cg D (M n ) and 
6"( 2 ^) 2 < C < b~( n+a f~ 2 ) 2 - The best constants are estimated by carefully studying the existence of 
positive solutions for the ODE (jig]) . 

Remark 2.14 Recently, Blanchet et al. in [4] studied a special case of inequality (j46| (a = b = 1, and 
a = 2) under the additional condition: 

{l + \x\ 2 f- x u{x)dx = Q, for [3 <^—^. (49) 



Note that we do not assume (|49p in Theorem 12. 13[ and that we have found the best constants for (3 < 0, a 
case that was left open in [3]. 

2.3 Improved Hardy and Caffarelli-Kohn-Nirenberg Inequalities 

In [9] Caffarelli-Kohn-Nirenberg established a set inequalities of the following form: 

(f R n \x\- bp \u\Pdx)i < C atb J Rn \x\- 2a \Vu\ 2 dx for all u G C^{R n ), (50) 

where for n > 3, 

-oo < a < a^, a < 6 < a + 1, and p = „_ 2+ 2 2 " (b _ o) ■ (51) 
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For the cases n = 2 and n = 1 the conditions are slightly different. For n = 2 

— oo < a <0, a < b < a + 1, and p = r^-, (52) 

and for n = 1 

— oo < a < o+5<6<o + l, and p = _ 1+2 2 ( . b _ a - ) . (53) 
Let -D*' 2 be the completion of C£°(R n ) for the inner product (it, u) = J*,„ |a;| _2a VM.Vwc?x and let 

S(a,6)= inf ^ ".' ' ; 2/p (54) 



ue-Dj' 2 \{0} (J an |a:|-ip|u|Pdx) 



denote the best embedding constant. We are concerned here with the "Hardy critical" case of the above 
inequalities, that is when b = a + 1. In this direction, Catrina and Wang [TT] showed that for n > 3 we have 
S(a, a+1) = ( "~ 2 2 a ~ 2 ) 2 and that S(a, a+1) is not achieved while 5(a, 6) is always achieved for a < 6 < a+1. 
For the case n = 2 they also showed that 5(a, a + 1) = a 2 , and that S(a, a+1) is not achieved, while for 
a < b < a + 1, S(a, 6) is again achieved. For n = 1, S(a, a + 1) = ( 1 + 2a ) 2 is also not achieved. 
In this section we give a necessary and sufficient condition for improvement of (|50[) with b = a+1 and n > 1. 
Our results cover also the critical case when a = which is not allowed by the methods of [S] . 

Theorem 2.15 Let W be a positive radial function on the ball B in R n (n > 1) with radius R and centered 
at zero. Assume a < . The following two statements are then equivalent: 

1. W is a Bessel potential on (0,R). 

2. There exists c > such that the following inequality holds for all u € C^{B) 

(H a , cW ) J B |a;|- 2a |Vw(a;)| 2 da; > (li^p 2 ) 2 j B \ x \-2a-2 u 2 dx + c j g | a! |-aa W -( a .j ti 2 da . j 

Moreover, ( — 2a ~ 2 ) 2 is the best constant and [3(W;R) — sup{c; {H a ^ c w)holds} , where f3(W;R) is the weight 
of the Bessel potential W on (0, R). 

On the other hand, there is no strictly positive W € C 1 (0,oo), such that the following inequality holds for 
all u £ Cg°{R n ), 

\x\- 2a \Vu(x)\ 2 dx > ( - - 2 ° - 2 ) 2 / \x\- 2a - 2 u 2 dx + c [ W{\x\)u 2 dx. (55) 



Proof: It suffices to use Theorems 12.11 and 12.91 with V(r) = r 2a to get that W is a Bessel function if and 
only if the pair (r~ 2a , W 0lC (r)) ^ s a Bessel pair on (0, R) for some c > 0, where 

W 0)C (r) = C~ 2 y 2 ) 2 r- 2 - 2a +cr- 2a W{r). 



For the last part, assume that (|55| holds for some W. Then it follows from Theorem 12.151 that for V = 
cr 2a W(r) the equation y"(r) + -tf + v(r)y = has a positive solution on (0, 00). From Lemma l2~7l we know 

that y is strictly decreasing on (0, +00). Hence, v y , £j > — ^ which yields y'(r) < ^, for some 6 > 0. Thus 
y(r) —* —00 as r — ► +00. This is a contradiction and the proof is complete. □ 

Remark 2.16 Theorem 12 . 151 characterizes the best constant only when £1 is a ball, while for general domain 
Q, it just gives a lower and upper bounds for the best constant corresponding to a given Bessel potential W. 
It is indeed clear that 

C BR (W)<C n (W)<C Bp (W), 

where Br is the smallest ball containing f2 and B p is the largest ball contained in it. If now W is a Bessel 
potential such that (3(W, R) is independent of R, then clearly (3(W, R) is also the best constant in inequality 
(Ha,cw) for any smooth bounded domain. This is clearly the case for the potentials Wu, P and Wu. P where 
P{W, R) = 4 for all R, while for W = 1 the best constant is still not known for general domains even for the 
simplest case a = 0. 
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Using the integral criteria for Bessel potentials, we can also deduce immediately the following. 

Corollary 2.17 Let Q be a bounded smooth domain in R n with n > 1, and let W be a non-negative function 
in C 1 (0, R —: sup^gao |x|] anda<^^. 

1. If \imini\n(r) sW{s)ds > — oo, then there exists a := a(fi) > such that an improved Hardy 
inequality (H. a ,w a ) holds for the scaled potential W a (x) := a 2 W / (a|x|). 

2. If lim ln(r) J sW(s)ds = — oo, then there are no a,c > 0, for which (H a ,w a c ) holds with W a ,c — 
cW{a\x\). 

By applying the above to various examples of Bessel potentials, we can now deduce several old and new 
inequalities. The first is an extension of a result established by Brezis and Vazquez [8] in the case where 
a = 0, and 6 = 0. 

Corollary 2.18 Let ft be a bounded smooth domain in R n with n > 1 and a < ^y 2 -- Then, for any b < 2a +2 

there exists c > such that for all u € Cq > (Q) 

/„ \x\- 2a \Vu\ 2 dx > (!Mp2)2 \ x \~2a-2 u 2 dx ( 56 ) 

Moreover, when £1 is a ball B of radius R the best constant c for which 156]) holds is equal to the weight 
/3{r 2a ' b ;R) of the Bessel potential W{r) = r 2a - b on (0, R]. 
In particular, 

J B \x\' 2a \Vu\ 2 dx > ("^2)2 J b \ x \-2a-2 u 2 dj , + ^ ^ \ x \-^ U 2 d X , 

(57) 

where the best constant Xb is equal to ZQUhJ n \Q,\~ 2 / n , where u> n and \Cl\ denote the volume of the unit ball 
and respectively, and z$ = 2.4048... is the first zero of the Bessel function Jq(z). 

Proof: It suffices to apply Theorem 12.151 with the function W(r) — r b+2a which is a Bessel potential when- 
ever b > —2a— 2 since then liminfo^o hi(r) J Q s 2a+1 W(s)ds > — oo. In the case where b = —2a and therefore 

W = 1, we use the fact that /3(1; R) = (established in the appendix) to deduce that the best constant is 

then equal to ZQUj' 2 / n \^\~ 2 ^ n ■ □ 

The following corollary is an extension of a recent result by Adimurthi et all [1] established in the case where 
a = 0, and of another result by Wang and Willem in [27 (Theorem 2) in the case k = 1. We also provide 
here the value of the best constant. 

Corollary 2.19 Let B be a bounded smooth domain in R n with n > 1 and a < . Then for every integer 

e((fc-l)-iime.) 

k, and p = (sup^g^ M)( ee )> we have for any u € Hq(Q,), 

J a \x\- 2 «\Vu\ 2 dx > (Mf^) 2 /„ j^dx + \ ZU Inj^iUU logVftr'dx. (58) 
Moreover, 4 is the best constant which is not attained in ffg(O). 

Proof: As seen in the appendix, Wk, P {r) — Y^=i ~pz ( 1°9 ^ jfy ) ^ &c is a Bessel potential on (0, R) 
where R = s\xp x&n \x\, and /3(Wk,p] R) — \- □ 

The very same reasoning leads to the following extension of a result established by Filippas and Tertikas [13] 
in the case where a = 0. 

Corollary 2.20 Let f2 be a bounded smooth domain in R n with n > 1 and a < ^i— . Then for every integer 
k, and any D > sup^gQ \x\, we have for u G H^(fl), 

In l j^dx > i^ 2 ) 2 J n j^dx + I ZZx L ^Xli}4)X 2 {^)...X^)\u\ 2 dx, (59) 

and j is the best constant which is not attained in Hq(Q). 
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The classical Hardy inequality is valid for dimensions n > 3. We now present optimal Hardy type inequalities 
for dimension two in bounded domains, as well as the corresponding best constants. 

Theorem 2.21 Let Q, be a smooth domain in R 2 and € f2. Then we have the following inequalities. 

• Let D > sup xen \x\, then for all u £ ifg(fi), 

J a \Vu\ 2 dx > \ YZx In j^xU l 4)xH l 4)---x?( l 4)H 2 dx (60) 

and i is the best constant. 

e((*-l)-time«) 

• Let p = (sup^gQ |a;|)(e e ), then for all u £ Hq(D.) 

/ n |V«| a «fa > IZUfn^illLlogU^dx, (61) 

and j is the best constant for all k > 1 . 

• If a < 2, then there exists c > such that for all u £ Hq(H), 

f n \Vu\ 2 dx>cf nJ ^dx, (62) 
and the best constant is larger or equal to (3{r a ; sup |x|). 

An immediate application of Theorem 1 2 . 1 1 coupled with Holder's inequality gives the following duality state- 
ment, which should be compared to inequalities dual to those of Sobolev's, recently obtained via the theory 
of mass transport [3J [TU] . 

Corollary 2.22 Suppose that fi is a smooth bounded domain containing in R n (n > 1) with R := 
sup^gQ M- Then, for any a < and < p < 2, we have the following dual inequalities: 

inf { jf \x\~ 2a \Vu\ 2 dx - ( W ~ 2 2 °~ 2 ) a J \x\- 2a - 2 \u\ 2 dx; u £ Cg°(n), = lj 

r f \x\~ 2a v i 
3 General Hardy-Rellich inequalities 

Let £ C R n be a smooth domain, and denote 

Co r (0) ={»£ Cq (fi) : v is radial and supp v cO}, 

H™(fl) = {u£ H™(Cl) : u is radial}. 
We start by considering a general inequality for radial functions. 

Theorem 3.1 Let V and W be positive radial C 1 -functions on a ball B\{0}, where B is a ball with radius R 
in W 1 (n > 1 ) and centered at zero. Assume J Q r n-iy/ r s dr — oo and lim r ^o r a V(r) = for some a < n — 2. 
Then the following statements are equivalent: 

1. (V, W) is a Bessel pair on (0, R). 

2. There exists c > such that the following inequality holds for all radial functions u £ C^° r (B) 
(RRv,cw) J B V(x)\Au\ 2 dx > cf B W(x)\Vu\ 2 dx + (n - 1) J B (X& - ^ffl)\Vu\ 2 dx. 
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Moreover, the best constant is given by 

P(V, W; R) = sup {c; (HRy jC vy) holds for radial functions} . (63) 
Proof: Assume u S C§° r (B) and observe that 



V(x)\Au\ 2 dx = niu n { V^u^r^dr + in-l) 1 V^-^r^dr + 2(n - 1) / V{r)uu r r n - 2 dr}. 
Jo Jo r Jo 

Setting v = u r , we then have 

V(x)\Au\ 2 dx = [ V(x)\Vu\ 2 dx+(n-l) f _ Y^Ml)\ v \* dx , 



Thus, (HRyjy) for radial functions is equivalent to 

V{x)\\Jv\ 2 dx > I W(x)v 2 dx. 

J B 

Letting x{r) = v{x) where \x\ = r, we then have 

V{r){x'{r)fr n - l dr> I W {r)x 2 {r)^- 1 dr. (64) 
Jo 

It therefore follows from Theorem 12. II that 1) and 2) are equivalent. □ 



By applying the above theorem to the Bessel pair 

V(x) = \x\- 2m and W m (x) = V(x) [( "~ 2 2 m ~ 2 ) 2 |x|- 2 + W{x)) 

where W is a Bessel potential, and by using Theorem 12. 9[ we get the following result in the case of radial 
functions. 

Corollary 3.2 Suppose n > 1 and m < ^y^. Let Br C R" be a ball of radius R > and centered at zero. 
Let W be a Bessel potential on (0, R). Then we have for all u € C^ r {BR) 

Moreover, ( " + 2 2m ) 2 and j3{W;R) are the best constants. 
3.1 The non-radial case 

The decomposition of a function into its spherical harmonics will be one of our tools to prove the corre- 
sponding result in the non-radial case. This idea has also been used in [21]. Any function u £ (rj) could 
be extended by zero outside il, and could therefore be considered as a function in C x '(R n ). By decomposing 
u into spherical harmonics we get 

u = ^T=o u k where u k = f k (\x\)(p k (x) 

and (ip k (x)) k are the orthonormal eigenfunctions of the Laplace-Beltrami operator with corresponding eigen- 
values Cfe ~ k(n + k — 2), k > 0. The functions fk belong to C^fi) and satisfy f k (r) = 0(r k ) and 
f'(r) — 0(r k ~ r ) as r — > 0. In particular, 

99 = 1 and f = nuJ ^-y J 9Br uds = ^- / N=1 u{rx)ds. (66) 

We also have for any k > 0, and any continuous real valued functions v and w on (0, oo), 

V(\x\)\Au k \ 2 dx = f V (\x\)(Af k (\x\)-c k ^^) 2 dx, (67) 



R" 



and 



f W(\x\)\Vu k \ 2 dx= [ W(\x\)\Vf k \ 2 dx + c k f W(\x\)\x\~ 2 fldx. (68) 

Jr™ J R n J R n 
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Theorem 3.3 Let V and W be positive radial C 1 -functions on a ball B\{0}, where B is a ball with radius 
R in W 1 (n > 1) and centered at zero. Assume J r n-lyt r \ d>T = 00 an d linir— >o r a V(r) = for some 
a < (n-2). // 

W ( r ) ~ + - V rr( r ) > f° r < r < i?, (69) 

then the following statements are equivalent. 

1. (V, W) is a Bessel pair with f3(V, W;R) > 1. 

2. The following inequality holds for all u G Cfi°(B), 

(HIW) j B V(x)\Au\ 2 dx > J B W(x)\Vu\ 2 dx + (n - 1) J B (^ - Y^f)\Vu\ 2 dx. 

Moreover, if [3(V,W; R) > 1, then the best constant is given by 

[3(V, W; R) = sup {c; (HRy, cly ) holds) . (70) 

Proof: That 2) implies 1) follows from Theorem 13. II and does not require condition (|69[) . To prove that 1) 
implies 2) assume that the equation {By,w) has a positive solution on (0, R\. We prove that the inequality 
(HRv,w) holds for all u 6 C^°(B) by frequently using that 

FWI^'Wpr™- 1 ^ > W(j")a: 2 (r)>" n ~ 1 *- for all x E C x (0,i2]. (71) 
Indeed, for all n > 1 and /c > we have 

V(x)\Au k \ 2 dx = — f V (x)(Af k (\x\)~c k ^l) 2 dx 



V{r)(f>'(r) + ^fUr) - J-^-fr^dr 



V{r){f^{r)fr n - X dr + (n - l) 2 / ^(r)(/[(r)) 2 r"" 3 dr 

Jo 

+<3 f R V(r)f k (r)r n ' 5 +2(n-l) [ V{r)f'^r)f k (r)r- 2 
Jo Jo 

-2c fc /"* V(r)/£'(r)/ fc (r)r n - 8 dr - 2c fc (n - 1) /" V(r)/£(r)/ fc (r)r n - 4 dr. 
Jo Jo 

Integrate by parts and use ([66]) for k — to get 

V(x)\Au k \ 2 dx = { R V{r){f k \r)) 2 r n - l dr + (n-l + 2 Ck ) ( R V(r)(f k {r)fr n ^dr (72) 
Jo Jo 

+ (2c fc (n-4) + c 2 ) [ R V(r)r n - 5 f k (r)dr-(n-l) [ V r (r)r n - 2 (f' k ) 2 (r)dr 
Jo Jo 

- c fc (n-5) [ R V r (r)f 2 (r)r n -Ur - c k f ' v rr (r)f 2 {r)r n ^dr. 



Now define <?fc(r) = and note that g k {r) — 0(r k 1 ) for all k > 1. We have 



^W(/fcW) 2 r™- 3 = / y(r)(^(r)) 2 r"- 1 dr+ / 2y(r) 5fe (rK(r)r"- 2 dr + / V{r) g 2 {r)r n ^dr 

Jo Jo 

'^(s'tWlV- 1 *-- (n-3) [ R V(r)g 2 (r)r n - 3 dr- f* V r {r)g 2 (r)r n - 2 dr 



18 



Thus, 

Vl.Hfl.infr" \ : I W(r)Jt(r)r"-\lr-[n-:V) I V (r) f k {r)r n ~ b dr - I K(r)/2(r)r"" 4 dr. (7:i) 



R r R r R r R 



Substituting 2c k L V( r )(f'k( r )) 2rn 3 in G3) by its lower estimate in the last inequality (|75|) , we get 



1 



pR pR 

V(x)\Au k \ 2 dx > / W(r){f' k (r)) 2 r n - 1 dr+ / W{r)(f k (r)) 2 r n - 3 dr 
nuj n J Rn Jo Jo 

rR rR 

+ (n-1)/ V(r)(f k {r)) 2 r n - 3 dr + c k (n~l) / F(r)(/ fc (r)) 2 r"- 5 dr 
Jo Jo 

- { R V r (r)r n - 2 {f k ) 2 (r)dr-c k {n-l) ^ V r {r)r n -\f k ) 2 {r)dr 

Jo Jo 

r R 

+ c k (c k - (n- 1)) / ^(r)r"- 5 /IWdr 
Jo 

+ c fc / V(0 + HMO _ Vrr (r))f k (r)r n - 3 dr. 



The proof is now complete since the last term is non- negative by condition Note also that because of 

this condition, the formula for the best constant requires that f3(V, W; R) > 1, since if W satisfies (|()9l then 
cW satisfies it for any c > 1. □ 

Remark 3.4 In order to apply the above theorem to the Bessel pair 

V(x) = \x\- 2m and W m , c (x) = V(x) [( "~ 2 ^"~ 2 ) 2 |a;|- 2 + cW(x)] 

where W is a Bessel potential, we see that even in the simplest case V = landW / mjC (x) = { ! ^ L ) 2 \x\^ 2 +W{x) 1 
condition (J5U) reduces to (■ 2l j^) 2 |^| _2 + W(x) > 2 1 cc | 2 , which is then guaranteed only if n > 5. 
More generally, if V(x) — \x\~ 2m , then in order to satisfy (|6"9")l we need to have 



-in + 4) - 2\/n 2 -n+1 -in + 4) + 2\/n 2 -n + 1 

< m < , (74) 



and in this case, we have for m < and any Bessel potential W on Br, that for all u £ Cq°(.Br) 



b r M 2m - V 2 ' J Br \x\^ JBr Fr - 

Moreover, (™i^™) 2 and 0(W;R) are the best constant. 

Therefore, inequality ([75)) in the case where m — and n > 5, already includes Theorem 1.5 in [24\ as a 
special case. It also extends Theorem 1.8 in [24 where it is established under the condition 

-in + 4) + 2\/n 2 -n+1 

< m < — (76 

6 

which is more restrictive than (|74|) . We shall see however that this inequality remains true without condition 
(|74")l . but with a constant that is sometimes different from ( " + 2 2m ) 2 in the cases where (174")) is not valid. For 
example, if m = 0, then the best constant is 3 in dimension 4 and || in dimension 3. 

We shall now give a few immediate applications of the above in the case where m = and n > 5. Actually 
the results are true in lower dimensions, and will be stated as such, but the proofs for n < 5 will require 
additional work and will be postponed to the next section. 
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Theorem 3.5 Assume W is a Bessel potential on Br C R n with n > 3 ; then for all u € C^°(Bf{) we have 

r |vd 2 f 
\Au\ 2 dx>C(n) / -yi^-dx + 0(W; R) / W(x)\Vu\ 2 dx, (77) 
Jb r pi Jb r 

where C(3) = ||, C(4) = 3 and C(n) = ^ /or a// n > 5. Moreover, C{n) and j3(W;R) are the best 
constants. 

In particular, the following holds for any smooth bounded domain fi in R n with R = sup,,, go \x\, and any 

ue H 2 (n)nH^(n). 



• For any a < 2, 



and for a = 0, 



\Au\ 2 dx>C(n) I Ll^Ldx + P(\x\ a ;R) / l -^-dx, (78) 



J n \Au\ 2 dx > C(n) J n ^fdx + %f n \Vu\ 2 dx, (79) 
the constants being optimal when fi is a ball. 



e(k-Um.ea) 

For any k > I, and p = R(e e ), we have 



(80) 



3 — 1 * — 1 

• For D > R, and Xi is defined as M06\) we have 

I \Au(x)\ 2 dx >_ C(„) I + i±i ^(M«(M,. ^(M^ (S1) 

Moreover, all constants appearing in the above two inequality are optimal. 
Theorem 3.6 Let W{x) — W{\x\) be radial Bessel potential on a ball B of radius R in R n with n > 4, 



id such that 



W r {r) _ A 



W 



Hardy-Rellich inequality holds: 



'jyy = — + f(r), where f(r) > and \im r ^o r f ( r ) =0. If A < n — 2, then the following 



, .., , n 2 {n-A) 2 f u 2 , ,n 2 (n A 2) 2 . . nN f W(x) ., , 
Au\ 2 dx> \ ' / dx + {— + ± L )(3(W;R) / -j^u 2 dx, (82) 

B iD J B FI 4 4 J B 



Proof: Use first Theorem 13 . 5 1 with the Bessel potential W, then Theorem 12. 151 with the Bessel pair 
x\- 2 , \x\- 2 (^^\x\- 2 + W), then Theorem El with the Bessel pair (W, ( "~ X ~ 2)2 )\x\~ 2 W) to obtain 



\Au\ 2 dx > C(n) [ dx + [3(W, R) [ W{x)\Vu\ 2 dx 

B JB \ x \ JB 



> C(n) ( " 4)2 / ^-dx + C{n)l3{W,R) [ ^j-u 2 + p(W, R) / W(x)\\7u\ 2 dx 

4 JB FI JB \ x \ 

4 J B M 4 is FI 

2 

Recall that C(n) = \- for n > 5, giving the claimed result in these dimensions. This is however not the case 

when n = 4, and therefore another proof will be given in the next section to cover these cases. 

The following is immediate from Theorem 13.51 and from the fact that A = 2 for the Bessel potential under 

consideration. 

Corollary 3.7 Let Q be a smooth bounded domain in WL n , n > 4 and R — sup^g^ |x|. Then the following 
holds for all ueH 2 (fi) n H 1 (fi) 
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1. If p — R(e e ) and log^ l \.) is defined as U05\) , then 

/ n iM^« LP-^^t L^k^r^ <83 » 

j — 1 i — 1 

2. If D > R and Xi is defined as M06\) , then 

Theorem 3.8 Let W\(x) and Wzfe) be two radial Bessel potentials on a ball B of radius R in R n with 
n > 4. If a < 1, then there exists c(a,i?) > such that for all u G H 2 (B) n Hq(B) 

n — 2a — 2 f u 2 I u 2 

+c( ) 2 J^^-^ dx + cf3(W 2 ;R) J^W 2 (x)j-^dx, 

Proof: Here again we shall give a proof when n > 5. The case n = 4 will be handled in the next section. 
We again first use Theorem 13. 51 (for n > 5) with the Bessel potential |a;| _2a where a < 1, then Theorem 12. 151 
with the Bessel pair (|:r|~ 2 , \x 

|-2((n_4L| 2 .|-2 + W y^ then again 

Theorem 12.151 with the Bessel pair 

(|a;|- 2a , \x\- 2a {{ n - 2 £- 2 ) 2 \x\- 2 + W) to obtain 

\Au\ 2 dx > ^ / ^^dx + P(\x\- 2a -R) [ P^-dx 



D 

> 



n 2 r 


|V«| : 


4 1b 


\x\ 2 


n 2 (n — 


4)2 
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n 2 (n — 


4)2 
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/ ^ dx + ^-PiW^R) f W^x^dx + ^x^-R) f ^dx 
Jb \x\ 4 Jb \x\ Jb \x\ 



> j BW d X+ - mi ,R) jW l{x ) wdx 

n — 2o — 2 r v 2 r v 2 

+(3(\x\- 2a ; R)( ) 2 j B j^dx + [3(\x\~ 2a - R)(3(W 2 ; R) W 2 {x)^alx. 

The following theorem will be established in full generality (i.e with V(r) = r~ m ) in the next section. 

Theorem 3.9 Let W{x) = W{\x\) be a radial Bessel potential on a smooth bounded domain fl in R™, n > 4. 
Then, 

f a \Au(x)\ 2 dx - J n ^L dx _ £ J q W {x)u 2 dx > #||u||^, 2(n) u G H 2 (n). 

3.2 The case of power potentials \x\ m 

The general Theorem 13.31 allowed us to deduce inequality (|85|) below for a restricted interval of powers m. 
We shall now prove that the same holds for all m < IL ^-- The following theorem improves considerably 
Theorem 1.7, Theorem 1.8, and Theorem 6.4 in [2"4"j . 

Theorem 3.10 Suppose n> 1 and m < , and let W be a Bessel potential on a ball Br C R n of radius 
R. Then for all u G C^(B R ) 

r |Vul 2 
alx + (3(W;R) / W^jl—i-da;, (85) 
Jb r \ x \ 




where 

[ r i^idx 

-;ueC™(B R )\{0} 

Moreover, (3(W; R) and Q>m,n are the best constants to be computed in the appendix. 
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Proof: Assuming the inequality 



^ &n,m f I in i rtdX) 



Br M 2m - n ' m J BR M 2m+2 ' 



holds for all u £ C x> (Bfi), we shall prove that it can be improved by any Bessel potential W . We will use 
the following inequality frequently in the proof which follows directly from Theorem 12.151 with n=l. 



R -i r R. r R 

„af f rf„\\2 . ^ , a 1 \2 / n-2 t 2 



r 



(f(r)) 2 dr>(^—) 2 r a - 2 f 2 (r)dr + P(W;R) r a W(r)f 2 (r)dr, a > 1, (86) 



o 2 j j 



for all / £ C°°(0,R), where both (^-) 2 and /3(W;R) are best constants. 

Decompose u £ C^(Br) into its spherical harmonics T,^ =0 Uk, where u k — fk(\x\)(fk(x). We evaluate 
J k = Jr" ~^$r dx in tnc following way 

h = [ R r n - 2m -\rt{r)) 2 dr + [(n - l)(2m + 1) + 2c fc ] r n ~ 2m -*{f' k ) 2 dr 
Jo Jo 

r R 

+c k [c k + (n-2m- 4) (2m + 2)] / r"- 2ro - 5 (/ fc (r)) 2 dr 

Jo 

> /3(W0 / fl r »-2„ l -i vt/(a;)(/ , )2(ir + [( !i+2m )2 + 2cfe] f R rn ^ 3{f r fdr 

Jo z Jo 

f R 

+c k [c k + (n - 2m - 4)(2m + 2)] / r"- 2m - 5 (/ fe (r)) 2 dr 

Jo 

> /3(W) / r^-'WixXfifdr + an^ [ r n ~ 2m - 3 \f k f 'dr 

Jo Jo 

+/W[( ! ^^) 2 + 2c fc - a n>m ] j\ n - 2m - 3 W(x)(f k ) 2 dr 



i C ^ 4 ) 2 [(^4r^) 2 + 2c fe - a„, m ] + c fe [ Cfe + (n - 2m - 4)(2m + 2)]) /"* r- 2 ™- 5 (/ fc (r)) 2 dr. 



o 



Now by (|115|) we have 



, n-2m-4 , n + 2m , 

(( ^ ) '( 2 + 2Cfe ~~ + Cfe + - 2m - 4)(2rn + 2)J > c k a n ^ m , 

for all fe > 0. Hence, we have 

h > a n ^ m [ R r n - 2m - 3 (fl) 2 dr + a n ^c k [ R r n - 2m - 5 (f k (r)) 2 dr 
Jo Jo 

+(3(W) / r n - 2m - 1 W(x)(f k ) 2 dr + (3(W)i(——) 2 + 2c k -a n , m } / r n ~ 2m - 3 W{x)(f k fdr 
Jo z Jo 

> a n , m [ R r n - 2m - 3 (f k ) 2 dr + a n:m c k f" r n - 2m ~ 5 (f k (r)) 2 dr 
Jo Jo 

+0(W) I* r n - 2m - l W(x)(f k ) 2 dr + [3{W)c k r n ~ 2m ~ 3 W(x)(f k ) 2 dr 
Jo Jo 
\Vu\ 2 j f Jlr . ,|Vw |2 



B R \x\ 2m+2 ^ ' J Br " '\x\ 2m 



dx + 0(W) / W{xy—^dx. 



Moreover, it is easy to see from Theorem I2.1l and the above calculation that (3(W; R) is the best constant. 
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Theorem 3.11 Let f2 be a smooth domain in R n with n > 1 and let V € C 2 (0, R =: sup^gn \ x \) be a 
non-negative function that satisfies the following conditions: 

V r (r) < and j" rn J v(r) dr = - rn _l VAr) dr = +oo. (87) 

There exists Ai, A2 £ R such that 

r -$$ + \i >0 on (0,R) and lim ^ + A x = 0, (88) 

r^l+A 2 >0 on (0,R) and lim + A 2 = 3 (89) 

and 

(|(n - Ai - 2) 2 + 3(n - 3)) 7(r) - (n - 5)rK0) - r 2 K r (r) > for all r e (0, R). (90) 
TTien i/ie following inequality holds: 

,m 1 \ 1 a l2 , .(n-Ai-2) 2 , ^(n-Ai-4) 2 /" 2 , 
^(|x|)|Au| 2 ^ > i L + („_i))V 1 L "V u 2 dx 



(n-l)(n-A 2 -2) 2 /" K-(|x|) 
4 7 n |x|3 



,2 



cte. (91) 
Proof: We have by Theorem 12.91 and condition (1901) . 

— f V(x)\Au k \ 2 dx = [ R V(r)(rt(r))\- n - 1 dr + (n-l + 2c k ) [ R V(r)(f' k (r)) 2 r n - 3 dr 
nun J R n J Jo 

+ (2c fc (n - 4) + ci) f R V(r)r n - 5 f k (r)dr - (n - 1) /"* K.(r)r"" 2 (/0 2 (r)dr 
Jo Jo 

rR rR 

- Cfe(n — 5) / V r {r)fl(r)r n ~ A dr - c k I V r r(r)f£(r)r n ~ 3 dr 
Jo Jo 

' V(r)(^'(r)) V^dr + (n - 1) T V(r)(./£(r)) 2 r n - 3 dr 



(»- I) l' R V r {r)r n - 2 (f k ) 2 {r)dr 
+ r, j ( ( ±( n - Ax - 2) 2 + 3(n - 3)) V(r) - (n - 5)rV r (r) - rV„(r)j A 2 (r)r"" 5 dr 



The rest of the proof follows from the above inequality combined with Theorem 12.91 □ 

Remark 3.12 Let V(r) = r~ 2m with m < ^j^-- Then in order to satisfy condition (|9H)) we must have 

— 1 — V 1+ ^" jj _ < m < 2=^. Under this assumption the inequality (f9Tj) gives the following weighted second 
order Rellich inequality: 

|Ad 2 ,(n + 2m)(n-4-2m), 9 f u 2 

dx > — A L Y / j-^-ridx. 



B \x\ 2m ~ 4 ' J B \x\ 2m + 4 



In the following theorem we will show that the constant appearing in the above inequality is optimal. More- 
over, we will see that if m < — 1 — V 1+ ^ n jj ^ then the best constant is strictly less than ( ("+ 2m )fo 4 ~ 21 "^ ) 2 , 



This shows that inequality (|9T|) is actually sharp. 
Theorem 3.13 Let m < and define 

Pn,m= inf r 5 - ■ (92) 

„eC»(B)\{0} J^^L^dx 

. (n + 2m) (n — 4 — 2m) , . , , (n + 2m)(n — 2m — 4) 
/9 n ,m = (- ^-j -) + mm {fc(n + fc-2)[fc(n + fc-2) + ^ ^- i]}. 

4 fc— 0,1,2,... z 

Consequently the values of nm are as follows. 



Then 
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_ ( {n + 2m) (n - 4 - 2m) ^ 2 



2. //f-3<m<-l- 1)2 , toen 



/3„, m = ( (^ + ^)(n-4-2m) )2 + (n _ 1)[(n _ ^ + (n + 2m)(n - 2m - 4) ] 

5. J/jfe := "- 2 2 Tt ~ 4 e N, then 

P n , m = ( (n + 2m)( 7 4 - 2m) ) 2 + *(» + * - 2)[* (n + fc - 2) + (" + 2m)(n-2 ro - 4) ] 
^ // fc < "- 2 2 m - 4 < fc + 1 /or some k G N, then 

where 

[ , w , , (n + 2m)(n — - 2m — 4) n , s/ lr/ , (n + 2m)(n — 2m 

a(m, n, k) = min lk(n + k- 2)[k(n + k-2)+ y —— ^ '-],(k + l)(n + k — l)[(fe + l)(n + fc — 1) + n 



2 2 

Proof: Decompose u e C oo (Bji) into spherical harmonics £^ Wfc, where Ufc = fk(\x\)(fk(x). we have 



-/ ^Sf da; = / r^-\f^r)fdr+[(n-l)(2 m +l) + 2c k } [ r n -^- 3 (f k fdr 

ljJ n J R n \X\ Jo JO 

+ Cfe [ Cfe + (n-2m-4)(2m + 2)] / r n - 2m - 5 (f k (r)) 2 dr 

Jo 



^ (n + 2m){n - 4 - 2m) ^ 



2 Jo 



by Hardy inequality. Hence, 



^/ ,x , (n + 2m)(n — 4 — 2m) . 9 . . , „. ri/ , „. (n + 2m)(n — 2m — 4) nl 

/?n,m > S(n, m, fc := ^ J -±- i 2 + min {fc n+fc-2 fc n+fc-2 +- \ }• 

4 fe=o,i,2,... 2 

To prove that (3 n ,m is the best constant, let fc be such that 

in + 2m)(n — 4 — 2m). 9 , . , „, r , , , „ N (n + 2m)(n — 2m — 4), . . 

/3„, m = i ii- ^) 2 + fc(n + fc - 2)[fc(n + fc - 2) + ^ ^- ^. (93) 



Set 



u=\x\-^ ± + m+ ^k(xM\x\), 



where ifk{x) is an eigenfunction corresponding to the eigenvalue Cfc and ip{r) is a smooth cutoff function, 
such that < tp < 1, with ip = 1 in [0, |]. We have 

= (- {H ± 2m)in A - 4 - 2m) - c k + 6(2 + 2m + £ )) 2 + 0(1). 



Let now e — > to obtain the result. Thus the inequality 



|Au|2 >/3 f -^d, 



b„ M 2m -""' m 7 Bi? M 2m+4 
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holds for all u£C$°(B R ). 

To calculate explicit values of (3 n ,m we need to find the minimum point of the function 

(n + 2m)(n-2m-4) > 
f(x)=x(x-\ ), x>0. 

Observe that 

^ (n + 2m)(7i-2m-4) ) = Q 

To find minimizer k G N we should solve the equation 

k> + (n-2)k+ in + 2m)( - n - 2m - A) =0. 

4 

The roots of the above equation are x\ = " + 2 2m and X2 = "~ 2 ™~ 4 . l) follows from Theorem 13.111 It is 

easy to see that if to < —1 — then Xl < q. Hence, for to < — 1 — * L the minimum of 

the function / is attained in x-i. Note that if m < — 1 — V 1+ ^" jj ^ then B(n,ml) < B(n,m,0). Therefore 
claims 2), 3), and 4) follow. □ 

The following theorem extends Theorem 1.6 of [24] in many ways. First, we do not assume that n > 5 or 
to > 0, as was assumed there. Moreover, inequality ([M]) below includes inequalities (1.17) and (1.22) of [53] 
as special cases. 



Theorem 3.14 Let m < ^-g— and let W(x) be a Bessel potential on a ball B of radius R in R n with radius 



n— 4 

R. Assume = — ~ + f(r). where f(r) > and lim r _,o r f( r ) = 0. Then the following inequality holds 

for all u £ C^{B) 

|Am| 2 ^ f u 2 j 

B \x\^ dX ~ ^ m J B \x\^ dX 

w (ti + 2to) 2 in -2m- A - 2) 2 N f W(x) , , N 

+P(W;R) C ' + [ ~ - A — ) / jzAh^dx. (94) 



nu) n J Rn \x\ 2m 



Proof: Again we will frequently use inequality ([86]) in the proof. Decomposing u S C yo (Bn) into spherical 
harmonics E^Lgitfc, where Uk = fk(\x\)ifk(x), we can write 

' t |AUfc|2 ^ = I"' r n - 2m - 1 (tt(r)fdr+[{n-l)(2m + l) + 2c k ] f r n - 2m -*{f' k ) 2 dr 
Jo Jo 

+c k [c k + (n - 2m- 4)(2m + 2)] / r n - 2m ~ 5 {f k (r) fdr 

Jo 

> C^^f [ R r n - 2m - 3 (f k ) 2 dr + (3(W,R) C r-^W{x){f' k fdr 
1 Jo Jo 

+c k [c k + 2{ n ~ X ~ A f + (n - 2m - 4) (2m + 2)] /* r n - 2m - 5 (f k (r)) 2 dr, 
* Jo 
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where we have used the fact that c k > to get the above inequality. We have 



J n J R™ \x\ 2m 



+P{W-R) {n + 2 A m)2 [ R r n - 2 ™- 3 W(x)(f k ) 2 dr 



4 

R 

+/3(W;R) I r n ^ 2m " 1 W(x)(f k ) 2 dr 

^R 







> Pn, m I r n ~ 2m -*(f k ydr 

+ P(W;R)(^±f^ + ("- 2 "-*- 2 ) a )jfV»"- *W( X )(f k fdr 
^HL f ^S^dx 



0(W;R) (n + 2m) 2 (n — 2m — A - 2) 2 f W(x) 2 
nu) n 4 4 J B \x\ 2m+z 

by Theorem 12.91 Hence, holds and the proof is complete. □ 

Theorem 3.15 v4sswme —1 < m < and Zei W{x) be a Bessel potential on a ball B of radius R and 
centered at zero in R n (n > 1). Then there exists C > such that the following holds for all u G C^°(B): 

\Au\ 2 (n + 2m) 2 (n-2m-4) 2 f u 2 

- dx > Ta / TZH^TZ dx ( 95 ) 



B \x\ 2m ~ 16 J B \x\ 2m + 4 

+p{W;R) {n+ l m? J B ^^u 2 dx + ( 3(\x\ 2m ;R)\H\ Hi . (96) 

Proof: Decomposing again u <E C q x '(Br) into its spherical harmonics S^ Mfe where u k — f k (\x\)ip k (x), we 
calculate 



J n JR 



2 rR rR 

-dx = / r n - 2m - 1 (f' k '{r)) 2 dr+[(n-l)(2m+l) + 2c k ] / r n - 2m - 3 (f' k ) 2 dr 
Jo Jo 

f R 

+ c fe [c fe + (n-2m-4)(2m + 2)] / r n - 2m - 5 (f k (r)) 2 dr 

Jo 

> ( !i ^) 2 j ( fl r«- 2 ™- 3 (/[) 2 dr + / 3(|.| 2 ™; J R)^r«- 1 (/D 2 dr 

R 

n-2m-3( fl \2 , 



+ c k / r n - 2m -%f' k Ydr 
Jo 

> (" + 2 ^) 2 (^ ~ 2m ~ 4 ) 2 ^ r n-2m- 5(/)c) 2 jr 

( n + 2w ) 2 [ R W (r)r n - 2m - 3 (f k ) 2 dr 



4 

+ /5(|x| 2m ;i?) I* r^iftfdr + c^x^-R) f ' r n - 3 (f k ) 2 dr 
Jo Jo 

(n + 2ro) 2 (n- 2m- 4) 2 /" u 2 . 



/?(^;i?).(n + 2m) 2 f 2j , fl .. |2m 



16nw„ H 2m + 4 
Hence (ESI) holds. □ 



We note that even for m = and n > 4, Theorem 13.151 improves considerably Theorem A. 2. in [2]. 
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4 Higher order Rellich inequalities 

In this section we will repeat the results obtained in the previous section to derive higher order Rellich 
inequalities with corresponding improvements. Let W be a Bessel potential, (5 n ^ m be defined as in Theorem 
EIH and 

,,(n + 2m) 2 (n-2m-A-2) 2 x 
<r n , m = 0(W; R)(± r ^ + ^ '-). 



For the sake of convenience we make the following convention: FJ a, = 1. 

i=l 

Theorem 4.1 Let Br be a ball of radius R and W be a Bessel potential on Br such that ^ = — + f(r), 
where f(r) > and lim r ^o r f( r ) — 0. Assume m G N, 1 < I < m, and 2k + 4m < n. Then the following 
inequality holds for all u € Cq°(Br) 

f \A m u\ 2 , ±1 a f |A m -' M | 2 , ^ a f Wix)^™-*- 1 ^ 2 , 

-j^T^ -11 /3 ™^+ 21 7^ | x |2fc+4i ^ + 2^ CT - fc + 2 * 1 1 A».*+^-2 | x |2fc+^+2 ( 9? ) 

Proof: Follows directly from theorem 13. 141 □ 

Theorem 4.2 Let Br be a ball of radius R and W be a Bessel potential on Br such that = —7 + 

where f(r) > and lim r ^o r f( r ) — 0. Assume m € N , 1 < I < m, and 2k + 4m + 2 < n. Then the following 
inequality holds for all u € Cq°(Br) 



f |VA m d 2 , ,n-2fc-2. 2 U: f |A m -'u| 2 , 



i=0 

'-1 '-1 r- ttt/ w \m— i— 1..I2 



+ I ^ ) 2_^On,k+2 l+ l [\_P n ,k+2j-l / 1 |2fc+4»+4 ^ 



i=0 j=l 
/" IA" l d 2 

+ /3(W;iJ) / W(x)L—L dx (98) 
^Sj, Fl 

Proof: Follows directly from Theorem 12.151 and the previous theorem. □ 
Remark 4.3 For k — Theorems \4-l\ and \4-^\ include Theorem 1.9 in [2\l as a special case. 

Theorem 4.4 Let Br be a ball of radius R and W be a Bessel potential on Br such that = —7 + 

where f(r) > and lim r ^o r f{ r ) = 0. Assume m € N, 1 < I < m — 1, and 2k + 4m < n. Then the following 
inequality holds for all u G Cq°(Br) 



.n-2k-A. 2 tX ^ y Ty(x)|A m - 4 - 2 u| 2 J 

+ On,fc( 2 ) 2^ a ^ k + 2l + 2 YlP n ' k + 2 n U|2fe+4i+6 dX 

„■ — n „- — 1 J Br I I 



1=0 j=i 

+ f3(W;R)a n . k / W(x Y 1 dx + 0{W;R) I W{x) 1 ' " da; (99) 

Jb b Fl JBj, Fl 

where a n ^ m is defined in Theorem VS.KA 

Proof: Follows directly from Theorem |3~TU1 and the previous theorem. □ 
The following improves Theorem 1.10 in [24] in many ways, since it is assumed there that I < ~ n + s + 2 ^ n ' 2 -n+i 



and 4m < n. Even for k — 0, Theorem 14.51 below shows that we can drop the first condition and replace the 
second one by 4m < n. 
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Theorem 4.5 Let Br be a ball of radius R and W be a Bessel potential on Br such that . Assume m € N, 
1 < I < ra, and 2k + 4m < n. Then the following inequality holds for all u G C^(Br) 

da; ^ 11 J / MZk+Il dx ( 10 °) 



u„ \A 2k tl 4 •//>• I s 



i=l i=l Js « 11 

' tl a , +2 . 2 (n-2k-Aj) 2 f , JA^ttP 
+ /3(^;fi)^a„ |fc+2i _ 2 n 4 



J_ 'Zl n (n—Ob — A^ 2 r |Y7A"i-i„|2 

'.Br 

l (n-2k-4j) 2 f Tr ,,_JA m -*u| ; 

i=i j=i 
where a n ^ m are the best constants in inequality i85\) . 

Proof: Follows directly from Theorem 13. 101 □ 

5 Appendix (A): The class of Bessel potentials 

The Bessel equation associated to a potential W 

(B w ) y" + lv' + W{r)y = 

is central to all results revolving around the inequalities of Hardy and Hardy-Rcllich type. We summarize 
in this appendix the various properties of these equations that were used throughout this paper. 

Definition 1 We say that a non-negative real valued C 1 -function is a Bessel potential on (0, R) if there 
exists c > such that the equation (B c w) has a positive solution on (0, R). 
The class of Bessel potentials on (0, R) will be denoted by B(0, R). 

Note that the change of variable z(s) = y(e~ s ) maps the equation y" + ^y' + W(r)y — into 
(B' w ) z" + e- 2s W{e~ s )z{s) = 0. 

(101) 

On the other hand, the change of variables ip(t) — e , — - maps it into the nonlinear equation 

y \ e ) 

(B 1 ^) ^'(t)+^ 2 (t)+e- 2t W(e- t ) = 0. 

(102) 

This will allow us to relate the existence of positive solutions of (Bw) to the non-oscillatory behaviour of 
equations (B' w ) and (B 1 ^). 

The theory of sub/supersolutions -applied to (B' w ) (See Wintner already yields, that if (Bw) 

has a positive solution on an interval (0,R) for some non-negative potential W > 0, then for any W such 
that < V < W, the equation (By) has also a positive solution on (0,R). This leads to the definition of 
the weight of a potential W <E 23(0, R) as: 

(3(W;R) = sup{c > 0; (B c w) has a positive solution on (0, R)}. (103) 

The following is now straightforward. 

Proposition 5.1 1) The class 23(0, R) is a closed convex and solid subset o/C 1 (0,i?). 
2) For every W € 23(0, i?) ; the equation 

(Bw) V" + },V' + f3(W;R)W(r)y = 

has a positive solution on (0,R). 

The following gives an integral criteria for Bessel potentials. 
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Proposition 5.2 Let W be a positive locally integrable function on R. 

1. If liminf ln(r) /T sW(s)ds > — oo, then for every R > 0, there exists a := a(i?) > smc/i i/iai </ie 
scaled function W a (x) := a 2 W(ax) is a Bessel potential on (0,R). 

2. If lim In(r) fJ sVK(s)ds = — oo, then there are no a,c > 0, for which W a c = cW(a\x\) is a Bessel 

i — >o 

potential on (0,R). 

Proof: This relies on well known results concerning the existence of non-oscillatory solutions (i.e., those 
z(s) such that z(s) > for s > sufficiently large) for the second order linear differential equations 

z"(s)+a(s)z(s)=0, (104) 

where a is a locally integrable function on R. For these equations, the following integral criteria are available. 
We refer to [TBI El I2H1 HS1 [3D] ) among others for proofs and related results. 

i) If lim sup t _ >00 t J t °° a(s)ds < j, then Eq. Q104p is non-oscillatory. 

ii) If liminft^oo t J t °° a(s)ds > j, then Eq. (|104|) is oscillatory. 

It follows that if liminf ln(r) J Q r sW(s)ds > —oo holds, then there exists S > such that (Bw) has a positive 

solution on (0, S). An easy scaling argument then shows that there exists a > such that W a (x) := a 2 W(ax) 
is a Bessel potential on (0, R). The rest of the proof is similar. □ 

We now exhibit a few explicit Bessel potentials and compute their weights. We use the following notation. 

logW(.)=log(.) and logW (.) = log(log^ k -^ (.)) for k > 2. (105) 

and 

Xi(i) = (l-log(i)r 1 , X fc (t)=Xi(X fc _i(t)) fc = 2,3,..., (106) 
Theorem 5.1 Explicit Bessel potentials 

1. W = is a Bessel potential on (0, R) for any R > 0. 

D. T/ie Bessel function J is a positive solution for equation (Bw) with W = 1, on (0, Zq), where z = 
2.4048... is the first zero of Jo. Moreover, z$ is larger than the first root of any other solution for (Bi). 
In other words, 

for every R > 0, 

/3(l;fl) = ||. (107) 

3. If a < 2, then there exists R a > such that W(r) — r~ a is a Bessel potential on (0, R a ). 

e(fc-«me») 

4- For each k > 1 and p > R(e e ), the equation (Bi Wk ) corresponding to the potential 

W k , p (r) = Ttj =l Uj where U 3 (r) = ^(^=1 log™ ^ 

has a positive solution on (0, R) that is explicitly given by <Pk,p(f) — (Ili=i log^-)?. On the other 
hand, the equation [B lw +xu ) corresponding to the potential \Wk, P + At/fc has no positive solution 

for any A > 0. In other words, Wk tP is a Bessel potential on (0, R) with 

0{W k ;p,R) = yoranyk>l. (108) 
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5. For each k > 1 and R > 0, the equation (Si^ r ) corresponding to the potential 

W KR {r) = E* =1 ^ where U 3 {r) = £X?(£)X 2 2 (£) . . . XjUf £)*?(£) 
ftas a positive solution on (0, i?) f/iaf is explicitly given by 

Mr) = (Xi(£)x a (~) ■ ..Xk.^x^yi. 

On the other hand, the equation (Bi^ t R +\U h ) corresponding to the potential jWk.R + XUk has no 
positive solution for any A > 0. In other words, Wk,R is a Bessel potential on (0, R) with 

I3{W KR] R) = \ for any k > 1. (109) 
Proof: 1) It is clear that cp(r) — —log(j^r) is a positive solution of (Bq) on (0, R) for any R > 0. 

z 2 

2) The best constant for which the equation y" + + cy = has a positive solution on (0, R) is where 
zo = 2.4048... is the first zero of Bessel function Jq(z). Indeed if a is the first root of the an arbitrary solution 
of the Bessel equation y" + + y(r) = 0, then we have a < z . To see this let x(t) = aJ (t) + bY (t), where 
Jo and Yq are the two standard linearly independent solutions of Bessel equation, and a and b are constants. 
Assume the first zero of x(t) is larger than zq. Since the first zero of Yq is smaller than zq, we have a > 0. 
Also b < 0, because Yo(t) — » — oo as t — > 0. Finally note that Yq(zo) > 0, so if b < 0, then x(zq + e) < for e 
sufficiently small. Therefore, 6 = which is a contradiction. 

3) follows directly from the integral criteria. 

4) That <pk is an explicit solution of the equation [Bx w ) is straightforward. Assume now that there exists 
a positive function ip such that 

<p'(r) + rip" (r) 



zEi(n^ w ?)" 2 +(i+^(n^ w s 



^ v ' j=l i=l i=l 



Define f(r) = > 0, and calculate, 



p'{r) + rp"(r) _ p' k (r) + rp'^r) , f'(r) + rf'(r) f(r) ' 



<p(r) p k (r) f(r) /W^nLV(f) 



E 



Thus, 



If now f'(a n ) = for some sequence {a™}^! that converges to zero, then there exists a sequence {(3 n }^i 
that also converges to zero, such that f"(/3 n ) = 0, and f'((3 n ) > 0. But this contradicts (|110[) . which means 
that / is eventually monotone for r small enough. We consider the two cases according to whether / is 
increasing or decreasing: 

Case I: Assume f(r) > for r > sufficiently small. Then we will have 

(rf(r))'^ 1 



< 



E 



Integrating once we get 

f{r) > r-^—. — , 

rnLiog j (£) 
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for some c > 0. Hence, lim r ^o f( r ) = — 00 which is a contradiction. 
Case II: Assume f'(r) < for r > sufficiently small. Then 

(r/'(r))'>v- 1 



E 



rf'(r) -ttrn; =1 log J (^)' 
Thus, 



f'(r) > — : , (111) 

for some c > and r > sufficiently small. On the other hand 



Since /'(r) < 0, there exists I such that f(r) > i > for r > sufficiently small. From the above inequality 
we then have 

&/'(&) - af (a) < -^K—, — j — ). 

From (|lllj) we have lim a ^o a/'( a ) = 0. Hence, 

V(&) < V 



n; =1 io g j (f) 

for every b > 0, and 

/'M < - A/ 



rnf=iW(^)' 



for r > sufficiently small. Therefore, 

lim /(r) = +00, 

r—tO 

and by choosing Z large enouph (e.g., I > ■§ ) we get to contradict (|111[) . 

The proof of 5) is similar and is left to the interested reader. □ 

6 Appendix (B): The evaluation of a n m 

Here W6 evctliicitc tti6 best constants d n m which appear in Theorem 13. 101 
Theorem 6.1 Suppose n > 1 and m < ^y^. Then for any R > 0, the constants 

a n . m = inf I J r Bn gg - ; u € C^(B R ) \ {0} 

are given by the following expressions. 
1. For n = 1 

• if m € (— 00, — |) U [— |, — 5], i/ien 

,l + 2m. 2 
oi,m = ( — g — ^ 



i/ — I < m < — I , i/iera 



r/ n + 2m, 2 ( ("-4-2m)(n+2m) +2)2 

ai, m = mm{( ) , ( rt _ 4 , 2m )2 | }• 
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2. Ifm = then 

a m ,n = min{(n - 2) 2 , n - 1}. 

3. J/n > 2 andm< -("+4)+2V^-n +i > ^ ^ = (2 ^2m)2_ 

^. If2<n<3 and -("+4)+2V^~ n+i < m < or n > 4 a7ld < m < «^2 ; fen 

/ (n-4-2m)(n+2m) + ?l — l) 2 

5. For n > 4 and -(n+4)+2Vn"-n +i < m < r^i define k * = ^ - |)(n - 2)]. 
• Ifk*< 1, then 



a 



^ (n-4-2m)(n+2m) + ?l — l) 2 



• For £;* > 1 i/ie interval (mj := ("+ 4 )+ 2 V" 2 n+ T^ m 2 ._ a_4) C(m fr e divided in 2k* — 1 subintervals. 
For 1 < k < k* define 



1 _ 2{n - 5) - A/(n - 2) 2 - 12fc(fc + ra - 2) 

TO fc : - c > 



, 2(n - 5) + V(n - 2) 2 - 12fc(fc + n - 2) 

mi := - . 

k 6 

Ifm G (mj,m}] U [m 2 ,mg)] 7 then 

^ (n-4-2m)(n+2m) + ?l — l) 2 



( »-4- 2m ) 2 + n - 1 
• For k > 1 and m G ( m L TO fe+i] U [ TO fe+i' m fe)' 

( (n-4-2m)(n+2m) + + fc _ 2 )) 2 ( (n-4-2m)(n+2m) + (fc + ^ + fe _ ^2 

a„, m - mm{ , }. 

For m € (m^, , m 2 . ), £/ien 

_ , (n-4-2m)(n+2m) + fc* ( n + fc* - 2)) 2 ( (n-4-2m)(»+2m) + (fc , + !)( n + fc * _ X ))2 

a n , m - mm{ , ( n -4-2m )2 + (A; . + 1)(n + fc . _ }• 

Proof: Letting V{r) = r~ 2m then, 

W (r) - ^P- + - V„(r) = ( C~ 2 T~ 2 ) 2 - 2 - 4m - 2m(2m + l)) r - 2m - 2 . 

r z r 2 



In order to satisfy condition we should have 



-(n + 4) + 2V« 2 -n+1 ^ ^ -{n + 4) + 2^n 2 - n + 1 . . 

< 771 < . (112) 

6 6 

So, by Theorem 13.31 under the above condition we have a 7lym = ( " + 2 2m ) 2 as in the radial case. 

For the rest of the proof we will use an argument similar to that of Theorem 6.4 in [24] who computed a n , m 

in the case where n > 5 and for certain intervals of m. 

Decomposing again u G C^(Br) into spherical harmonics; u = T,^ =0 Uk, where Uk = fk(\x\)tpk{x), one has 
T^dx = / N- 2 ™(/n|x|)) 2 dx+((n-l)(2m+l) + 2 Cfe ) / \x\~ 2m - \f' k fdx (113) 
+ c fe (c fe + (n-4-2m)(2m + 2)) / \x\- 2m -\f k ) 2 dx, 
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lV '" t] \dx= I |*r 2 — \f k )*dx + n I ;,| '-"" '{h^b: (INI 



i„ lxl 2m + 2 

One can then prove as in [24] that 

a«,m = niin {A(k, m, n); k S N} (115) 

where 



( („ -4-2„ : ) ( , + 2m) + Cfc)2 _ 
( „-4-2 m )2+Cfc 

and 



4(fc,m,n) = L__^ 3 __±£SL i f m= n^ (n6) 



A(k, m, n) := c fc if m = and n + A: > 2. (117) 

Note that when m — and n + k > 2, then c& ^ 0. Actually, this also holds for n + k < 2, in which case 
one deduces that if to = then 

2 n + 2m 2 



a n ,m = min{(n - 2) = ( ) , (n - 1) = c x } 

which is statement 2). 

The rest of the proof consists of computing the infimum especially in the cases not considered in [Mj . For 
that we consider the function 

i (n-4-2m)(n+2m) 

/(*)=' 



( n-4-2m ) 2+a; 



2 

It is easy to check that /'(x) = at x\ and #2, where 

(n — 4 — 2m) (n + 2m) 
xi = (118) 

(n-4-2m)(-n + 6m + 8) 
a; 2 = . (119) 

Observe that for for n > 2, < ^j^. Hence, for m < both xi and X2 are negative and hence 
a n<m = (Hi^™)2. Also note that 



(n + 4) - 2Vn 2 -n + 1 ^ n-8 , „ ^ 1 
< /or a/i n > 1. 



6 

Hence, under the condition in 3) we have a n , TO = ( " + 2 2? " ) 2 . 

Also for n = 1 if m < — | both critical points are negative and we have ai )TO < ( 1+2t " ) 2 , Comparing 
j4(0, m, n) and m, n) we see that A(l, to, n) > A(0, m, n) if and only if (|112p holds. 
For n = 1 and — | < to < — g both and 22 are positive. Consider the equations 

(2m + 3)(2m + l) 
x{x -l) = x x = - , 

and 

(2m + 3)(6m + 7) 

x[x -l)=x 2 = • 

4 

By simple calculations we can see that all four solutions of the above two equations are less that two. Since, 
A(l,m, 1) < A(0, m, 1) for to < —2, we have ai, m < min{^4(l, to, 1), A(2, to, 1)} and 1) follows. 
For n > 2 and 2^ < m < we have x\ > and X2 < 0. Consider the equation 

(n - 4 - 2m)(n + 2m) 
x(x + n - 2) = xi = -- -. 



Then 2m + 4 n an d — ( 2m +") are solutions of the above equation and both are less than one. Since, for n > 4 



n-2 -(ra + 4) + 2 v / n 2 - n+1 
— — > — ■ 



and A(l,m, n) < A(0, m, n) for m > ("+ 4 )+2Vn n+\ ^ ^ e j-,^ constant is equal to what 4) claims. 
5) follows from an argument similar to that of Theorem 6.4 in |24j . 
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